The small-scale dynamo: Breaking universality at high Mach numbers 
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The small-scale dynamo may play a substantial role in magnetizing the Universe under a large 
range of conditions, including subsonic turbulence at low Mach numbers, highly supersonic turbu- 
lence at high Mach numbers and a large range of magnetic Prandtl numbers Pm, i.e. the ratio 
of kinetic viscosity to magnetic resistivity. Low Mach numbers may in particular lead to the well- 
known, incompressible Kolmogorov turbulence, while for high Mach numbers, we are in the highly 
compressible regime, thus close to Burgers turbulence. In this study, we explore whether in this 
large range of conditions, a universal behavior can be expected. Our starting point are previous 
investigations in the kinematic regime. Here, analytic studies based on the Kazantsev model have 
shown that the behavior of the dynamo depends significantly on Pm and the type of turbulence, 
and numerical simulations indicate a strong dependence of the growth rate on the Mach number of 
the flow. Once the magnetic field saturates on the current amplification scale, backreactions occur 
and the growth is shifted to the next-larger scale. We employ a Fokker-Planck model to calculate 
the magnetic field amplification during the non-linear regime, and find a resulting power-law growth 
that depends on the type of turbulence invoked. For Kolmogorov turbulence, we confirm previous 
results suggesting a linear growth of magnetic energy. For more general turbulent spectra, where the 
turbulent velocity Vt scales with the characteristic length scale as ui ex. £^ , we find that the magnetic 
energy grows as (t/Ted) 2 ^^ 1- ^- 1 , with t the time-coordinate and T e d the eddy-turnover time on the 
forcing scale of turbulence. For Burgers turbulence, ■& = 1/2, a quadratic rather than linear growth 
may thus be expected, as the spectral energy increases from smaller to larger scales more rapidly. 
The quadratic growth is due to the initially smaller growth rates obtained for Burgers turbulence, 
and thus implies longer timescales until saturation is reached. Similarly, we show that the char- 
acteristic length scale of the magnetic field grows as tVC-— *) - m ^ e general case, implying t 3//2 for 
Kolmogorov and t 2 for Burgers turbulence. Overall, we find that high Mach numbers, as typically 
associated with steep spectra of turbulence, may break the previously postulated universality, and 
introduce a dependence on the environment also in the non-linear regime. 



I. INTRODUCTION 

The small-scale dynamo has been suggested to operate 
under a large range of different conditions, including the 
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solar surface [T|, 0] , galaxies and galaxy clusters , the 
intergalactic medium |9| and the formation of the first 
stars and galaxies [10t - Tl7| . It thus operates on a large 
range of different conditions, concerning for instance the 
magnetic Prandtl number Pm, i.e. the ratio of kinematic 
viscosity v to magnetic resistivity rj, the Mach number of 
the turbulence M, i.e. the ratio of turbulent velocities 
to the sound speed, and, most likely related to the Mach 
number, the expected type of turbulence in the system. 

Most studies of the small-scale dynamo performed so 
far have focused on incompressible Kolmogorov turbu- 
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lence [l8|, assuming a scaling relation ui oc I 1 / 3 between 
turbulent velocity and length scale I. For Kolmogorov 
turbulence, it was previously concluded that the mag- 
netic energy grows exponentially in the kinematic regime 
[e.g. Il9l - t22j and linearly once the backreactions from the 
magnetic field become important (e.g. [23T425T ] ) . The 
latter was interpreted by Beresnyak [25J as evidence for 
universality of the small-scale dynamo, suggesting that a 
fixed fraction of the global turbulence dissipation rate is 
converted into magnetic energy. 

However, observations of turbulence in molecular 
clouds [e.g. Hi| H3] and numerical simulations of super- 
sonic turbulence [e.g. [28l - l3l| often reveal steeper turbu- 
lent spectra, typically inbetween the incompressible Kol- 
mogorov turbulence and the highly compressible Burgers 
turbulence [1H. So far, only a small amount of studies 
have investigated the turbulent dynamo in this regime. 
For instance, Haugen et al. [H[ provided the first study 
exploring the dependence of the dynamo on the Mach 
number in simulations of driven turbulence, and Balsara 
et al. [34l . l35l | , Balsara and Kim [36[ explored the ampli- 
fication of magnetic fields in turbulence produced from 
strong supernova shocks. The first systematic study cov- 
ering turbulent Mach numbers from 0.02 to 20 and differ- 
ent types of turbulence driving has been pursued by Fed- 
errath et al. [!t| , while the effect of a large range of dif- 
ferent turbulence spectra has been explored by Schober 
et al. [23| based on the Kazantsev model [l9j . 

We note that the small-scale dynamo has also been 
studied in the context of so-called shell models [38l - l4l| . 
The latter originate from shell models of hydrodynam- 
ical turbulence, which originally considered turbulence 
in 2D [42TO |. but were extended to 3D once a descrip- 
tion of kinetic helicity was obtained [45] . The first 2D 
MHD shell model has been derived by Frik 46], while 
3D models have been developed by Brandenburg et al. 
[13, Basu et al. [3, Frick and Sokoloff (H . M ore sophis- 
ticated processes such as non-local interactions [13, HI| , 
anisotropies [H[ and the Hall effect [H| have been incor- 
porated in more recent studies. These approaches allow 
to study both the evolution of the power spectrum as well 
as the saturated regime, and are highly complementary 
to the methods presented here. 

In the following, we will consider the small-scale dy- 
namo in the kinematic and non-linear regime, and present 
evidence from existing and new calculations suggesting 
a strong dependence on the magnetic Prandtl number, 
as well as the Mach number of the flow. In section 2, 
we summarize the evidence and indications for a non- 
universal behavior in the kinematic regime, which has 
been derived in previous studies. In section 3, we present 
the first exploration concerning different types of turbu- 
lence during the non-linear phase of the dynamo, where 
the backreaction of the magnetic field becomes impor- 
tant. We show that a linear growth is only obtained in 
the case of Kolmogorov turbulence, while steeper power 
laws result from turbulent spectra with > 1/3. We 
discuss the physical implications in section 4, and sum- 



marize our main results in section 5. 



II. NON-UNIVERSALITY IN THE KINEMATIC 
REGIME 

The small-scale dynamo is well-studied in the kine- 
matic regime, where an exponential growth of the mag- 
netic field is expected on the viscous scale. The growth 
rate of the magnetic field can be calculated in the 
framework of the Kazantsev model, assuming homoge- 
neous turbulence that is ^-correlated in time, or with 
3-dimcnsional magneto-hydrodynamical simulations. In 
this section, we discuss hints and evidence for a non- 
universal behavior in the kinematic regime. 

A. Indications for non-universality in the 
Kazantsev model 

The amplification of magnetic fields is governed by the 
induction equation, which is given as 

d t B = V x v x B + rjAB. (1) 

We assume in the following (B) = 0, although we note 
that scenarios considering (B) ^ have been recently 
explored by Boldyrev et al. [5J] , Malyshkin and Boldyrev 
[55L HH, H3 • I n the Kazantsev model, the velocity field 
and the magnetic field are decomposed into a mean field, 
denoted with brackets (), and a fluctuating component 
denoted with 5: 

v = (v) + 5v, B = (B)+ SB. (2) 

A central input is the correlation function of the turbu- 
lent velocity, which is <5-correlated in time and (in the 
absence of helicity) can be decomposed as 

(dv i (f 1 ,t)Svj(f 2 ,s)) = Tij(r)S(t - s), (3) 

T ij (r)=(5 ij -^)T N (r) + ^T L (r), 

with r = |ri — and T/v, Tl are the transverse and lon- 
gitudinal parts of the correlation function, respectively 
[58|] . The same definitions can be applied to the magnetic 
field, yielding a two-point correlation function My(r, t) 
with transverse and longitudinal components Mjv (r, t) 
and Ml(t, t). Unlike the velocity field, the magnetic field 
is always divergence-free, leading to the additional con- 
straint 



As the Kazantsev model assumes that the flow is 6- 
correlated in time, concepts such as viscosity or the mag- 
netic Prandtl number cannot be directly incorporated 
into the flow, as the turbulent velocity field is destroyed 
and regenerated at each instant, leaving no time for vis- 
cosity to act. However, it can be indirectly included by 
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adopting turbulent velocity spectra that become steeper 
below a given viscous scale t v . This is the approach em- 
ployed here. For a given relation of type 



(5) 



in the inertial range, the longitudinal correlation function 
of turbulence can be parametrized as [22j 



VL 
3 



r L (r)=^(i_( £ )^ 







< r < 4 
l v <r < L 
L <r, 



(6) 

with t v the viscous scale, L the driving scale of turbu- 
lence, V the turbulent velocity on scale L, Re = VLjv 
the Reynolds number of the gas and v the kinetic viscos- 
ity. Similarly, we have 



VL 
3 



T N (r) = ^(i_^)(r 




0+1 



< r < 4 
l v <r < L 
L <r, 

(7) 

with 0(0) = (21 — 380)/5. As we expect an exponential 
growth of the magnetic energy as a function of time, we 
make the following ansatz for the kinematic regime: 



M L (r,t) = 



1 



2Tt 



(8) 



Inserting © in the induction equation ([T]), one obtains 
the Kazantsev equation, which is of the same form as the 
quantum-mechanical Schrodingcr equation: 



Kdiff(r) 



dV(r) 



d 2 ? 



+ U{r)ijj{r) = -Tip(r). (9) 



In this framework, the amplification depends on the ef- 
fective potential U(r) in Eq. ©, which depends on the 
properties of turbulence via 



TJ( r \ = K diff _ ( K diff) 2 , ^ K dif 

1 >~ 2 4 Kdiff + r2 
Kdiff =»7 + T £ (0) -T L {r). 



2(t l - r N ) 



(10) 



As recently shown by Schober et al. [22}, this form of 
the potential also accounts for the effect of compressibil- 
ity by keeping terms related to V-w during the derivation. 
The equation can be solved using the WKB approxima- 
tion in the limit of Pm — > oo [19l - [22l |. For Kolmogorov 
turbulence, one obtains 



rV, Pm »i = 1.028 -Re 



1/2 



(11) 



In a recent study, analytical solutions based on the WKB 
approximation have been derived in the limit Pm <C 1 by 
Schober et al. [59j]. For Kolmogorov, they yield 



Jf,Pm«l 



0.0268^Rm 1/2 , 
Li 



(12) 



with Rm = VL/rj the magnetic Reynolds number, and r\ 
the magnetic diffusivity. We thus observe a fundamen- 
tal difference between the limiting cases Pm <C 1 and 
Pm 3> 1 in the kinematic regime: For Pm 3> 1, magnetic 
field amplification occurs predominantly on the viscous 
scale, corresponding to the most negative range of the 
potential. For Pm <C 1, on the other hand, the resis- 
tive scale becomes larger than the visous scale. Amplifi- 
cation on the viscous scale is thus not possible, and the 
strongest contribution is close to the resistive scale due to 
the short eddy-times. Correspondingly, the growth rate 
of the magnetic field depends on the Reynolds number 
Re for Pm 3> 1, and on the magnetic Reynolds number 
Rm for Pm -c 1 [see also H3] ■ 

The results can be generalized further for different 
types of turbulence. In the limit Pm ^> 1, one obtains 



J.Pm» 1 



(163 - 3041) (!-*)/(!+*) 
60 L 



(13) 



In the regime Pm <C 1, one finds a similar relation [59j, 



L 



(14) 



with the prefactor a defined through the quantities 

a(0) = 0(56 - 1030), (15) 
6(0) = 0(79 - 1570), (16) 



25 + J lib aid) + (6(0) - 25) 2 - 6(0) 
c(0) = ¥ . . — . (17) 



as 



«(0) 



C (0)' 



exp 



a(0) 



3a(0) 



7T (0-1) -2 . (18) 



A numerical evaluation shows that these coefficients are 
smaller by about two orders of magnitude in the limit 
Pm <C 1, assuming the same type of turbulence. This can 
be expected, as the amplification then occurs on larger 
scales, with larger eddy-turnover times. 

Similarly, also the type of turbulence reflected in the 
parameter may change the amplification rate by about 
an order of magnitude, in case of the same value of 
Rm. The most efficient amplification rate occurs for Kol- 
mogorov turbulence, = 1/3, while it is less efficient 
for highly compressible Burgers turbulence, = 1/2, for 
which the turbulent velocities decrease more rapidly with 
length scale. 

The Kazantsev model thus indicates that the behav- 
ior of the dynamo depends both on the type of turbu- 
lence and the magnetic Prandtl number. A potential 
restriction of the underlying model is the assumption 
of ^-correlated turbulence, although the characteristic 
timescales are certainly small compared to the dynam- 
ical time. To investigate the resulting uncertainties, we 
refer the reader to Schekochihin and Kulsrud 61 1. The 



main results from these considerations are thus the fol- 
lowing: 
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• The behavior of the small-scale dynamo depends 
sensitively on the value of Pm, and in particular 
whether Pm < 1 or Pm ^> 1. We note that there 
is a continuous transition at Pm ~ 1, as detailed 
by Bovino et al. [62j . 

• The adopted type of turbulence has a significant 
influence on the efficiency of magnetic field ampli- 
fication, as turbulent spectra with ?9 > 1/3 cor- 
respond to larger eddy-turnover times and smaller 
amplification rates. 



B. Results from numerical simulations 

Due to the numerical viscosity and resistivity, it is dif- 
ficult to perform magneto-hydrodynamical simulations 
with Pm significantly different from 1. However, a lim- 
ited range of Pm has nevertheless been explored. For in- 
stance, Haugen et al. [63| investigated magnetic Prandtl 
numbers between 0.1 and 30. For Pm < 1, they report 
that the miminum magnetic Reynolds number required 
for dynamo action, Rm c , scales as 



Rm, - 357rPm~ 1/2 . 



(19) 



We note that the factor ir in the above is due to their def- 
inition of the magnetic Reynolds number. They further 
report differences in the obtained power spectra, indicat- 
ing a steeper decrease on small scales for small values of 
Pm. 

Schckochihin et al. 64 1 and Iskakov et al. [65j report 
numerical simulations exploring the small-scale dynamo 
from Pm~ 0.017 up to Pm =1. In this regime, they 
find that even for constant values of Rm, the growth rate 
decreases with decreasing Pm. In particular, for Pm ~ 1 
and Rm ~ 830, they report a normalized growth rate of 
1.8, which decreases to 0.9 for Pm ~ 0.2 and the same 
magnetic Reynolds number. The simulations further 
indicate that the value of Rm c settles to a constant limit 
for Pm < 1, Re> 1 and Rm 3> 1, even though this case 
is hard to numerically explore. The scaling of the growth 
rate on Rm, on the other hand, has not been conclusively 
explored. 

All in all, simulations thus show that the growth rate 
depends on the magnetic Prandtl number even in the 
range Pm < 1. Another quantity which was shown to 
influence the dynamo is the Mach number of the gas. 
Haugen et al. [33| explored Mach numbers in the range 
of 0.1 — 2.1 and reported a clear dependence of the criti- 
cal magnetic Reynolds number for dynamo action on the 
Mach number Ai. For Pm ~ 5, they report Rm c ~ 25tt 
for M. < 1 and a rapid increase to Rm c ~ 45ir for M. > 1 . 
A similar behavior was found for Pm ~ 1, with critical 
values of ~ 407r and ~ 807r, respectively. 

A larger series of simulations has been reported by Fed- 
errath et al. (37j , exploring Mach numbers from 0.02 up to 
20, with compressive and solcnoidal forcing, respectively. 
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FIG. 1. A fit to the growth rates obtained by Federrath et al. 
[3?| for compressive and solenoidal forcing as a function of 
Mach number. These simulations correspond to Pm ~ 2. 



They performed a fit to the growth rate and saturation 
levels as a function of Mach number, using the function 



f(M)= po 



M P1 



P2 



M? 3 + Pi 



p 5 )M r 



(20) 



The fit coefficients for the different cases are given in 
Table U and the normalized growth rates are given in 
Fig. [U In the subsonic regime, their results indicate that 
the growth rate (normalized by the eddy-turnover time 
T e d on the forcing scale) strongly decreases with decreas- 
ing Mach number for compressive driving, while it is al- 
most constant at solcnoidal driving. At Ai > 1, there 
is an initial drop due to the appearance of shocks, but 
increases as Ai 1 ^ 3 at larger values. A similar dependence 
is reported on the saturation level, which is particularly 
high for solenoidal driving, and decreases in the regime 
of large, supersonic Mach numbers. 
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Fcomp[T cd ] 


(E m / Ek)sol 
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PO 


-18.71 


2.251 


0.020 


0.037 


pi 


0.051 


0.119 


2.340 


1.982 


v> 


-1.059 


-0.802 


23.33 


-0.027 


P3 


2.921 


25.53 


2.340 


3.601 


P4 


1.350 


1.686 


1 


0.395 


PB 


0.313 


0.139 





0.003 


Pa 


1/3 


1/3 









TABLE I. Fit coefficients reported by Federrath et al. [33] . 

We thus summarize the results from numerical simu- 
lations as follows: 

• The critical magnetic Reynolds number for dynamo 
action as well as the resulting spectra for the mag- 
netic field depend on the magnetic Prandtl number. 
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• Both the growth rates and the saturation levels of 
the dynamo depend significantly on the turbulent 
Mach number and the type of forcing that is em- 
ployed. 



III. NON-UNIVERSALITY IN THE 
NON-LINEAR REGIME 

The exponential growth phase will come to an end 
when the tension force of the magnetic field, B ■ VB, be- 
comes comparable to the inertial term of the flow, u- Vm. 
At this point, magnetic field amplification will stop on 
the scales that fulfill this condition, and continue to pro- 
ceed on larger scales. As discussed by Schckochihin et al. 
1231 , this condition translates to 



Br 



(21) 



where £ a denotes the smallest scale where amplification 
still occurs. In this regime, a linear growth of the mag- 
netic energy has been reported in previous studies, based 
on the assumption of Kolmogorov turbulence [e.g. [23l - l25| . 
In the following, we will generalize these investigations 
by employing a simplified toy model as well as a more 
sophisticated Fokker-Planck model previously suggested 
by Schekochihin et al. [HI. As a result, we will show 
that different types of power-law growth can be expected 
depending on the adopted type of turbulence. 

We further point out that in the non-linear regime, 
we expect the magnetic Prandtl number to play a less 
critical role, as the amplification scale of the magnetic 
field is now expected to be larger than both the viscous 
and the resistive scale, such that no strong dependence 
on Re or Rm can be expected. 

We note that the models considered in this section have 
previously been motivated in the context of the incom- 
pressible induction equation, given as 



d t B + v ■ \7B = B ■ W + r]AB. 



(22) 



However, they can be naturally extended into the com- 
pressible regime with the replacement 



B 



B 
P ' 



(23) 



Inserting this replacement as well as the continuity equa- 
tion, 



p = -V • (pv) , 



(24) 



it is straightforward to show that one obtains the com- 
pressible form of the induction equation, 

d t B + v-VB = B ■ W-B(V- v) + 7]AB, (25) 

equivalent to Eq. JT]). As long as the mean density (p) 
in the box is constant, a significant growth of the quan- 
tity (B I p) nevertheless implies a corresponding growth 



of the magnetic energy, assuming that the density distri- 
bution function will not change significantly over time. 
In the case of well-developed driven turbulence, one in- 
deed expects a characteristic log-normal density prob- 
ability distribution function, which naturally complies 
with these requirements j66M68| . Strictly speaking, the 
following considerations apply to the quantity B = B/p 
and W — W/p 2 , with W the magnetic energy. In the 
following, the ~ is however dropped for simplicity. 



A. First considerations based on a toy model 

In the toy model previously proposed by Schckochihin 
et al. [23|, the dominant fraction of the magnetic energy 
resides on the scale £ a , the smallest scale where magnetic 
field amplification still occurs (thus yielding the shortest 
amplification timescale). On that scale, the magnetic en- 
ergy is expected to be already close to saturation. The 
magnetic energy W(t) can thus be related to the ampli- 
fication scale l a by the approximate relation 



1 



W(t) ~ =(p)ui a(t) . 



(26) 



The magnetic energy is evaluated here at the mean den- 
sity (p) of the turbulent box, as we are interested only in 
the magnetic field amplification by shear. Adopting the 
eddy-turnover rate on the scale l a as the growth rate for 
the magnetic field, i.e. 



r(f) 



u i*{t) 



the magnetic energy evolves as 
d 



<lt 



W = T(t)W(t) - 2 v k 2 rms W(t) 



with 



and 



2 1 [°° 



dkk 2 M{t, k) 



M{t,k) = - I dQ~ k (\B(t,k)\^ 



(27) 



(28) 



(29) 



(30) 



Now, we have T(t)W(t) ~ {p)u\, t) /l a (t) =: e(t). Insert- 
ing in Eq. (J2HJ) yields 

jW = X e(t)-2 V k 2 rms (t)W(t), (31) 

where x is a constant of order unity. For Kolmogorov tur- 
bulence, the quantity e(t) = (p)uf / t \/£a(t) is a constant 
[l8j . In this case, and as long as magnetic energy dis- 
sipation is negligible, dW/dt = const, implying a phase 
of linear growth. In this limit, we obtain the result of 
Beresnyak [25|, where a constant fraction of the turbu- 
lence dissipation rate is converted into magnetic energy. 



6 



In the case of Burgers 
"■ 5 . In this case, the 



In the general case with ug a oc e(t) is however not 
constant, but varies as t^f~ l . 
turbulence, we thus obtain e oc 
growth of the magnetic energy is no longer linear, as the 
turbulent energy dissipation rate is not independent of 
scale! 



For comparison, 

3/2 



note that the quantity e 



pe S Q S /£, with esGS the specific energy density of subgrid- 
scale turbulence, is practically independent of t. It how- 
ever has a weak dependence on the Mach number, and 
a strong dependence on the type of forcing [69[. As the 
density fluctuations will however not contribute to the 
shearing, we will adopt e as the quantity of interest here. 

To quantify the expected behavior, wc need to solve 
Eq. (|26[) for l a . For this purpose, we recall that ui a is 
related to the turbulence driving scale L and the velocity 
V on that scale via 



Ug 



V 



From we thus obtain 



£a = L 



\(p)V 2 J 



1/(2*) 



(32) 



(33) 



We can now evaluate (|27| and (|28|) . yielding 



-W~ W 



( 2W 



1/(2*)' 



0-1 



x 



W /l+(*-l)/(2*)_ 



\(p)v 2 , 

(34) 

For Kolmogorov turbulence (0 = 1/3), we confirm that 
dW/dt = const, while in the more general case, this 
quantity will increase with increasing W. This can be 
intuitively understood, as the steep spectra for d > 1/3 
imply a more modest increase of the eddy-timescale with 
length scale, suggesting that the amplification rate re- 
mains larger when increasing the scale. We re-assess 
these results with the Fokker-Planck model below and 
explore the physical implications in more detail. 



B. Implications of the Fokker-Planck model for 
universality 

The starting point for our investigations is the Fokker- 
Planck model of Schekochihin et al. [23j]. Here, the time- 
evolution of the magnetic-energy spectrum is given as 



d t M 



d_ 

dk 



D(k) 



dM 
~dk 



V(k)M + 2T(t)M ~2rjk 2 M, 

(35) 

with the diffusion coefficient D(k) = T(t)k 2 /5 and the 
drift velocity in k-space V(k) = AT(t)k/h. We recall 
that the magnetic-energy spectrum M is related to the 
magnetic energy W via 



W(t) 



dkM(t,k). 



(36) 



To describe the evolution in the nonlinear regime, 
Schekochihin et al. [23| postulated the following expres- 



r(t) = ci 



MO 



dkk 2 E{k) 



1/2 



W(t) = c 2 / dkE(k). 

Jk B (t) 



(37) 
(38) 



The constants c\ and c-i arc of order unity, E(k) is the 
hydrodynamic energy spectrum neglecting the influence 
of the magnetic field, and the wave vector k s (t) is de- 
fined via Eq. (|38[) . It corresponds to the smallest scale 
where amplification efficiently occurs. As input for the 
Fokker-Planck model, we require an energy spectrum of 
the turbulence. As before, we assume that the velocity 
in the incrtial range scales as 



ut<x£*. 



(39) 



The hydrodynamic energy spectrum is then approxi- 
mately given as 



E{k) 



Ct€ 2/3 k -20- 





1 for k £ [kf, k v ] 
elsewhere, 



(40) 



with Ct a constant which depends on the type of tur- 
bulence, kf and k v the wave vectors describing the in- 
jection scale of turbulence and the viscous scale, respec- 
tively. The value of k v is set to enforce the condition 
e = 2v J^dkk 2 E(k). Unlikely in © and 0, we do 
not explicitly model the turbulent spectra in the viscous 
regime, as these no longer contribute during the non- 
linear stage. With these input data, Eq. (f38| can be 
evaluated as 



w(t) = [k 2 * 



We further introduce the quantities 

20 



W = c 2 dkE(k) = 
Jo 



-201 



L-20 _ , -20 



C2C t e 2 /* 

2d v ' 



(41) 

,(42) 
(43) 



We note that in the above expressions, the integral J °° dk 
corresponds to an integration from kf to k v , as the tur- 
bulent energy is non-zero only in this regime (see l40|) . 

Using these definitions, the wave vectors k v , k s and kf 
can be expressed as 



2d 



2i) 



-1/(20) 



[W{t) + W„]- 1/{2i> \ (44) 



/ 9,9 \ - 1 /( 2 ^ 



— - 



2d 



c 2 C t eV 3 



-1/(20) 



(46) 
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Integrating Eq. (|37|) now yields the following: 

• Ct£ 2/ 3 X . .-I V2 



r(t) = ci 



20 



,2-2,9 



(47) 



Substituting Eqs. (|4"4" l) -(|4l) ]l into (gTJ yields the expression 
'C t e 2 /3\ 1/2 /c 2 C t e 2 / 3 \^ 



r(t) = d 



2-2?? 



21? 



(48) 



(W(t) + Wj/) 1- ' - (1K + IK,) 1 "* 



1/2 



.(49) 



Considering turbulence models between Kolmogorov and 
Burgers, we have 1/3 < •& < 1/2. We further assume 
that W(t) <C Wo, implying that the magnetic field is far 
from saturation on the current amplification scale. In 
this case, we can neglect the second term in the square 
brackets. As we focus here on the non-linear regime, we 
can further neglect W v compared to W(t), and obtain 
the expression 



r(t) = ci 



g t e 2/3 
2-2?? 



1/2 



( c 2 C t e 



2/3 



V 2ti 



w W-l)/(?&)Q.y 
(50) 

As in our toy model, the growth of the magnetic energy 
thus scales as 

^-WocWWmocW^-WV. (51) 
at 

For Kolmogorov turbulence, the growth is thus linear, 
while it grows faster than linear for $ > 1/3. Integrating 
Eq. ([ST]) , we obtain 



(52) 



with 



C 



C t e 2 / 3 
2-2?? 



1/2 



C-2 



20 



;i-^)/(2i?) 



i 

20 



(53) 

From this expression, we already see that the energy 
grows linearly in t for Kolmogorov, while it grows as t 2 
for Burgers turbulence. For a physical interpretation, the 
normalization in terms of the eddy-turnover time T e d on 
the forcing scale is still required, which we perform below. 



IV. PHYSICAL IMPLICATIONS 

To explore the physical implications of the above- 
mentioned results, we now perform a normalization in 
terms of the eddy-turnover time T e d on the forcing scale 
kj 1 . For this purpose, we note that the expression within 
the central brackets of Eq. ([5B")) is identical to W v k v , and 
it is straightforward to show that 



(54) 




0.0001 0.0010 0.0100 

Time t [T ed ~ 



0.1000 



1.0000 



FIG. 2. The power-law growth of magnetic energy for differ- 
ent types of turbulence in the non-linear regime, following the 
evolution Eq. J56| for Re = 10 4 . 



Model and reference 


r? 


W oc 


la OC 


Kolmogorov 1181 


1/3 


t 1 


£3/2 


Intermittency of Kolmogorov turbulence [70] 


0.35 


t 1 - 077 


£1.54 


Driven supersonic MHD turbulence [28] 


0.37 


t 1.17 


£1.59 


Observation in molecular clouds [26] 


0.38 


£1.23 


£1.61 


Solenoidal forcing of turbulence [31] 


0.43 


t 1.51 


£1.75 


Compressive forcing of turbulence [311 


0.47 




£1.89 


Observation in molecular clouds [27] 


0.47 


£1.77 


£1.89 


Burgers turbulence [32] 


1/2 


t 2 


t 2 



TABLE II. The power-law behavior of the small-scale dynamo 
for different types of turbulence in the non-linear regime. 



If we normalize Eq. (|52"j) in terms of T e d ~ (kf^/Wo) 1 , 
we thus obtain 



(55) 



C 



C t e 2 / 3 
2-2?? 



1/2 



1 

2^ 



x Re M ( l - M )/( l +^)fc 2 *" 2 ^ /(l - ,?) W 2l ' _1 ' /(l_,?) 

Adopting a system of units with Wo = 1 and kf = 1, 
it is evident that Ef ~ 1, v(kf) ~ 1 and thus T e d ~ 1. 
From Eq. (|40|) . we also expect e ~ 1. In these units, our 
evolution equations simplifies as 

2i3/(l-i5) 

- I , (56) 



W(t) 



Wn 



c = 



t 
i 



2 - 2$ 



1/2 



1 

2i) 



Re 2tf(l-2tf)/(l-#) < 



We illustrate the behavior for the different types of tur- 
bulence in Fig. [5] for Re = 10 4 , and summarize the 
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power-law behavior in Table [TTJ The solution suggests 
that turbulence spectra closer to Kolmogorov saturate 
earlier (in terms of the eddy-turnover time on the forc- 
ing scale kf), and initially start at a higher value. The 
latter is fully consistent with our expectations for the 
kinematic regime, where the growth rates are higher for 
Kolmogorov turbulence, and a larger amount of mag- 
netic energy may build up before the non-linear regime is 
reached (due to the increased amount of turbulent energy 
that is available on the same scale). We note that in the 
final stage close to saturation, the evolution may start to 
deviate from the power-law behavior reported here, pro- 
viding a transition to the regime where W(t) = const. 

From the relation derived above, we further calcu- 
late the characteristic scaling of the current amplification 
scale l s as a function of time t. Adopting Eq. (ffij)) , we 
have W(t) ~ (p)uj, t) oc 0f, thus 



E(k)~k 



-2-6-1 



4(t) ocWVP^f) 0^1/(1- 



(57) 



For Kolmogorov turbulence, the characteristic length 
scale of the magnetic field thus grows as t 3 / 2 , while it 
grows as t 2 for Burgers turbulence. The results are sum- 
marized for all types of turbulence in Table (|TT|) . 

The power-laws derived here depend on the type of tur- 
bulence due to the different eddy-turnover timescales as a 
function of scale, as we sketch in Fig. ©. We summarize 
the main ingredients based on the toy model developed 
in section IIII Al 

Considering a driving scale L with a turbulence veloc- 
ity V on that scale, the ratio of the eddy-turnover times 
on scale I <C L for Kolmogorov and Burgers turbulence 



is given as 



t_K 



-1/3 



(£/L) 



1-1/2 



1/6 



(58) 



During the growth of the magnetic energy, the relevant 
length scale however shifts to larger scales. According 
to Eq. (|55|) . the ratio of the eddy timescales approaches 
unity for £ — > L. For Burgers turbulence, the magnetic 
field amplification is thus initially delayed with respect 
to Kolmogorov, and catches up later, resulting into the 
non-linear behavior and the power-law growth described 
here. 

Due to these results, it is clear that the growth rate of 
the dynamo is not a fixed fraction of the global turbu- 
lence dissipation rate, as previously proposed by Beres- 
nyak [25|. Due to the dependence on the turbulent spec- 
trum, such a consideration may only hold locally, i.e. on 
a given scale, where the growth rate of the field is in- 
deed related to the local eddy timescale. From a more 
global perspective, however, the turbulence dissipation 
rate changes as a function of scale for models different 
from Kolmogorov, such that the previously postulated 
universal behavior cannot be expected. From Eq. (f56|). 
it is further evident that the evolution depends on the 
Reynolds number of the gas, and that larger Reynolds 
numbers imply stronger magnetic fields at earlier times. 



Kolmogorov 




Burgers N N 

-2 



cf. Eq. (36) 



FIG. 3. A sketch of Kolmogorov vs Burgers turbulence. While 
the turbulent energy is considerably smaller for Burgers spec- 
tra (i? = 1/2) on small scales, it approaches the values for 
Kolmogorov turbulence (i? = 1/3) on larger scales. As a re- 
sult, the magnetic energy grows faster than linear for Burgers 
turbulence, as the growth rates gradually approach the Kol- 
mogorov values at later times. 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have explored both the kinematic 
regime of the small-scale dynamo, where an exponential 
growth of the magnetic energy is generally observed, and 
the non-linear regime, where backreactions start occuring 
on small scales and shift the amplification scale of the 
magnetic field to larger scales. 

In the kinematic regime, analytical studies based on 
the Kazantscv model suggest a fundamental dependence 
on the magnetic Prandtl number. In particular, for 
Pm <C 1, the growth rate of the dynamo is a function of 
the magnetic Reynolds number Rm, while for Pm ^> 1, 
it depends on the kinematic Reynolds number Re. In 
addition, the amplification rates significantly depend on 
the adopted type of turbulence. For Pm ^> 1, it scales as 
Re 1 / 2 for Kolmogorov turbulence and as Re 1//3 for Burg- 
ers turbulence. The same scaling relations, with a dif- 
ferent normalization, were found for Pm -C 1, with the 
replacement Re— >Rm. 

Numerical simulations confirm the dependence on Pm 
also in the range Pm ~ 1, and find a strong dependence 
of the growth rate and the saturation level on the turbu- 
lent Mach number M. and the type of turbulence forc- 
ing. Magnetic field amplification is particularly efficient 
for solenoidal forcing and low Mach numbers, but also oc- 
curs for high Mach numbers and solenoidal / compressive 
forcing. If the Mach numbers are very small, compressive 
forcing is hardly able to trigger magnetic field amplifica- 
tion, as the presence of density gradients are required for 
the production of solenoidal turbulence in this case. 
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To investigate the non-linear regime of the dynamo, 
we employed the Fokker-Planck model of Schekochihin 
et al. [23| and explored the effect of different turbulent 
spectra on the magnetic field amplification rate. We find 
that the previously known linear growth only occurs for 
Kolmogorov turbulence, while in the general case with 
ui oc i , we expect the magnetic energy to scale as 
i 2 */ C 1 -*). The energy growth is thus faster than linear, 
and may even become quadratic for Burgers turbulence 
({) = 1/2). However, we note that the growth rate is 
initially smaller for Burgers turbulence, as the turbu- 
lent energy available for amplification is initially much 
smaller on small scales. While magnetic field amplifica- 
tion is shifted to larger scales, the difference in the turbu- 
lent energy decreases, implying the reported power-law 
behavior as a function of time. 

We have further shown that also the scaling of the 
characteristic length scale £ a for magnetic field amplifi- 
cation depends on the turbulent slope. Specifically, we 
find a scaling as i 1 '^ - *', corresponding to t 3 / 2 for Kol- 
mogorov and t 2 for Burgers turbulence. The change of 
length scales proceeds thus in a fashion analogous to the 
inverse-cascade in case of helicity [e.g. [Tl|, [72J. The evo- 
lution of this quantity may thus provide another relevant 
diagnostic for a comparison with numerical simulations. 

Due to the above considerations, we point out that 
the non-linear stage of the small-scale dynamo does not 
generally correspond to converting a fixed fraction of the 
turbulence dissipation rate into mag netic energy, as pre- 
viously suggested by Beresnyak [231. While their results 
agree with our model for the case of Kolmogorov tur- 
bulence (low Mach numbers), steeper power laws may 



occur in the highly compressible regime. Universality in 
the sense of a uniform behavior under all conditions can 
thus not be expected. Nevertheless, we note that there 
are still universal laws governing the behavior of the dy- 
namo, which relate the growth of the magnetic energy 
to the eddy-turnover time on the current amplification 
scale. This quantity in general does depend on the Mach 
number and the type of turbulence involved, such that 
the breaking of universality is a result of the properties 
of different environments. We propose to explore such 
effects in further detail with numerical simulations to im- 
prove our understanding of such non-universal behavior. 
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